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ABSTRACT 


Paint Lyman-G (LyG) emitters become increasingly rarer towards the re-ionisati on epoch (z ^ 6 — 7). However, 
observations from a very large 5 deg^) Lya narrow-band survey at 2 ; = 6.6 ( Matthee et al. [2015 ) show that 
this is not the case for the most luminous emitters, capable of ionising their own local bubbles. Here we 
present follow-up observations of the two most luminous Lya candidates in the COSMOS field: ‘MASOSA 
and ‘CR7’. We used X-SHOOTER, SINPONI and PORS2 on the VET, and DEIMOS on Keck, to confirm 
both candidates beyond any doubt. We find redshifts of ^ = 6.541 and 2 : = 6.604 for ‘MASOSA’ and ‘CR7’, 
respectively. MASOSA has a strong detection in Lya with a line width of 386 ± 30kms“^ (PWHM) and 


with very high EWq (> 200A), but undetected in the continuum, implying very low stellar mass and a likely 
young, metal-poor stellar population. ‘CR7’, with an observed Lya luminosity of ergs“^ is the 

most luminous Lya emitter ever found at 2 ; > 6 and is spatially extended 16kpc). ‘CR7’ reveals a narrow 
LyG line with 266 ± 15kms“^ PWHM, being detected in the NIR (rest-frame UV; (3 = —2.3 ± 0.1) and in 


IRAC/Spitzer. We detect a narrow Heiil640A emission line (6cr, PWHM = 130±30 km s in CR7 which can 
explain the clear excess seen in the J band photometry (EWq ^ 80 A). We find no other emission lines from the 
UV to the NIR in our X-SHOOTER spectra (Heii/Oiii]1663A > 3 and Heii/Ciii]1908A > 2.5). We conclude 
that CR7 is best explained by a combination of a PopIII-like population which dominates the rest-frame UV and 
the nebular emission and a more normal stellar population which presumably dominates the mass. i75'7yWPC3 
observations show that the light is indeed spatially separated between a very blue component, coincident with 
LyG and Hell emission, and two red components 5 kpc away), which dominate the mass. Our findings are 
consistent with theoretical predictions of a PopIII wave, with PopIII star formation migrating away from the 
original sites of star formation. 


Subject headings: galaxies: evolution, cosmology: early Universe, cosmology: dark ages, re-ionisation, first 
stars. 


1. introduction 

The study of the most distant sources such as galaxies, 
quasars and gamma ray bursts offers unique constraints on 
early galaxy and structure formation. Such observations are 
particularly important to te st and refine models of galaxy for¬ 
mation and evolution (e.g. [Vogelsberger et al.]|2014 |Schaye| 
et al.||2015|) and to study the epoch of re-iomsation (e.g. 


Shapiro et al.|| 19941 [Sokasian et al.||2004 |Eurlanetto et al. 


2004HMcQuinn et al.||2006HIliev et al.||2066| ). Over the last 

two decades, considerable effort has been dedicated towards 
finding the most distant sources. More recently, and particu 
larly due to the upgraded capabilities of HST, multiple can 
didate galaxies u p to 2 ; ^ 8 — 11 (e.g. |Bou wens et al.|2011 
Ellis et al.|2013 ) have been found with deep broad-band pho 
tometry. However, spectroscopic confirmation is st ill limited 
to a handful of galaxies and quasars at 2 : > 6.5 (e.g. |Mortlock 


let al.||20lT| |Qno et aL]|2Q12t |Schenker et al.|[^14t |Einkei 
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Istein et aL||2Q13[ |Pentericci et al.||MT^ |Qesch et ar||2015| ), 

for both physical (galaxies becoming increasingly fainter) and 
observational reasons (the need for deep near-infrared expo¬ 
sures). At these redshifts {z > 6.5), the Lya line is virtually 
the only line available to confirm sources with current instru¬ 
ments. However, Lya is easily attenuated by dust and neu¬ 
tral hydrogen in the inter-stellar and inter-galactic medium. 
Indeed, spectroscopic follow-up of UV-select ed galaxies in¬ 
dicate that LyG is supp ressed at 2 ; > 7 (e.g. jCaruana et al. 
20141 |Tilvi et ah |2014| ) and not a single 2 ; > 8 Lya emit 


ter candidate has been confirmed yet (e. g. Sobral et al.|2009t 


Eaisst et ^|2Q14[ |Matthee et al. 2014| ). If the suppression 
of LyG is mostly caused by the increase of neutral hydrogen 
fraction towards higher redshifts, it is clear that 2 : ^ 6.5 (just 
over 0.8 Gyrs after the Big Bang) is a crucial period, because 
re-ionisation sh ould be close to complete at that redshift (e.g. 
|Ean et al.|2006| ). 


and confirming Lya emitters at 2 : 

^3 — 7 (e.g. ICowie & Hu 

1998 

IMalhotra & Rhoads|2004[ 

lye et al. |2006[ 

jMurayama 

et al. 

120071 |Hu et al.||2010l |Ouchi et al.||20101). 

The results 


show that the Lya luminosity function is constant from 2 ; ^ 3 
to 2 ; ^ 6, but there are claims that the number density drops 
from 2 : ^ 6 to 2 ; ^ 6.6 (e.g. Ouchi et^pOlOj Kashikawa 


et al.|201 1]) and that it drops at an even faster rate up to 2 : ^ 7 
(e.g. fshibuya et al.|2012 JKonno et al.|2014| ). Moreover, the 
^ ^ /function 


fact that the res t-frame UV luminosit y lunction declines from 
2 ; ^ 3 — 6 (e.g. Bouwens et al.|20T4] ) while the Lya luminos- 
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ity functio n (LF) is roughly c onstant over the same redshift 
range (e.g. [Quchi et~ar 2008b implies that the cosmic aver¬ 
age Lya e s^cape fraction is likely increasing, from ~ 5% at 
z ^ 2 (e.g. [Hayes et al.|20T0} Ciardullo et al.|2014]), to likely 
^ 20 — 30% around z ^ 6 (e.g. jCassata et al.||2015j ). Sur¬ 
prisingly, it then seems to fall sharply with increasing redshift 
beyond 2 : ^ 6.5. Current results could be a consequence of re¬ 
ionisation not being completed at 2 ; ^ 6 — 7, particularly when 
taken together with the decline in the fraction of L yman break 
selected galaxies with high EW Lya emission (e .g. jTilvi et al.j 
[20141 [Caruana et al.|2014l [Pentericci et al.|2014b . However, it 


is becoming clear that re-ionisation by itself is not enough to 
explain the rapid de cline of the fraction of strong Lya em itters 
towards 2 ; ^ 7 (e.g. [Dijkstra|2014[[Mesinger et al.|2015] ). 

It i s likely that re-iomsatio n was very heterogeneous/patchy 
(e.g. Pentericci et al.[[MT4) , with the early high density re¬ 
gions rMomsmgTirstTT^nowed by the rest of the Universe. 
If that were the case, this process could have a distinguish¬ 
able effect on the evolution of the hya luminosity function, 
and it may be that the luminous end of the luminosity func¬ 
tion evolves differently from the fainter end, as luminous hya 
emitters should in principle be capable of ionising their sur¬ 
roundings and thus are easier to o bserve. This is exactly what 
is found by [Matthee et al.[([2Q15[), in agreement with spectro¬ 
scopic results from Qno et al.[ ( |2012| . 


In addition to using hya emitters to study re-ionisation, 
they are also useful for identifying the most extreme, metal 
poor and young galaxies. Studies of hya emitters at 2 ; > 
2 — 3 show that, on average, these s ources are indeed very 
metal poor ([Finkelstein et ar][2Qllt [Nakajima et al.[[MT2 


Guaita et al.[[2013' ), prese nting high ionisation p arameters 
(high [Olll]/H/i/ line ratios; [N^ajima et al.[[2Q13) and very 
low typical dust extinctions (kg. [Qno etal.[20 


Given these 


observations, hya searches should also be able to find metal- 
free, PopIII stellar populations (since galaxies dominated 


by PopIII emit large amounts of hya photons, e.g. [Schae^ 
2002[[20()^ . However, so far, although some ca ndidates for 


PopIIIstelfar populations have been found (e.g. [Jimenez & 


Haiman[[2006t [Dijkstra & Wyithe[[2007| [Nagao k al.[[200S^ 


Kashikawa et al.|2012[[Cassata et al.|2013[), and some meta i 


poor galaxies have been confirmed (e.g. Prescott et al.|2009] ), 

they are all significantly more metal rich than the expected 
PopIII stars, and show e.g. Cm] and Civ emission. For 
example, when there is no evidence for the presence of an 
AGN and no metal lines, the short-lived Helll640A emis¬ 
sion line (the ‘smoking gun’ for PopIII stars in extremely high 
EW Lyman-a emitters without any metal emission line) was 
never detected with hig h enough EW (e.g. Nagao et al.|2008 ; 
[Kashikawa et al.|2Q12[ ). 

Until recently, hya studies at the epoch of re-ionisation 
have been restricted to the more numerous, relatively faint 
sources of Li^ya ^ IQ^ ^'^ergs"^ (with some exceptions, 
5.7 follow -up: Westra et al. [2006} [Lidman et al. 


e.^ _ ___ ___ 

[2()12[ and ‘Himiko’: [Quctii et al.[[2009 ). However, with the 
wide-field capabili ties of current instrum ents (including Hy¬ 
per Suprime-Cam; [Miyazaki et al.|2012[ ), the identification of 
luminous hya emitters will become increasingly easier. Re¬ 
cently, significant progress wa s made towards find ing lumi¬ 
nous hya emitters at 2 ; = 6.6 ( [Matthee et al.|2015] ), through 
a ~ 5 deg^ narrow-band survey, which resulted in the iden¬ 
tification of t he most luminous hy a emitters at the epoch of 
re-ionisation. [Matthee et al.[([2Q15|) re produced the hya lumi¬ 


nosity function of Quchi et al. ( 2010[ ) for relatively faint hya 


emitters at 2 ; = 6.6 for the UDS field, who find a decrease 
in th eir number density co mpared to lower redshifts. How¬ 
ever, [Matthee et al.[ ( [2015[ ) find that the luminous end of the 
2 ; = 6.6 LE resembles the 2 ; = 3 — 5.7 luminosity function, 
and is thus consistent with no evolution at the bright end since 
2 : ^ 3. Extremely luminous hya emitters at ^ ^ 6.6 are thus 
found to be much more common than expected, with space 
densities of 1. 5to:? X 10 ^ Mpc The results may mean 
that, because such bright sources can be observed at 2 ; ^ 6.6, 
we are witnessing preferential re-ionisation happening around 
the most luminous sources first. Such luminous sources may 
already be free (in their immediate surroundings) of a sig¬ 
nificant amount of neutral hydrogen, thus making their hya 
emission observable. Eurthermore, these sources open a new 
window towards exploring the stellar populations of the most 
luminous hya emitters at the epoch of re-ionisation even be¬ 
fore the James Webb Space Telescope (JWST) and ^ 30—40 m 
class telescopes (Extremely Large Telescopes, ELTs) become 
operational, as these are bright enough to be studied in un¬ 
precedented detail with e.g. HST, ALMA, VLT, Keck. 

Here we present spectroscopy of the two most luminous 
hya emitters found so far at the epoch of re-ionisation {z ^ 
7). This paper is organised in the following way. §2 presents 
the observations, the hya emitter sample, and the data reduc¬ 
tion. §3 outlines the details of the optical and near-infrared 
spectroscopic observations and measurements with the VLT 
and Keck data. §4 presents the discovery of the most lu¬ 
minous hya emitters and comparison with previous studies. 
§5 discusses SED fitting, model assumptions and how HST 
high spatial resolution data corroborates our best interpreta¬ 
tion of the data. §6 presents the discussion of the results. Ei- 
nally, §7 outlines the conclusions. A Hq = 70kms“^ Mpc“^, 
VLm = 0.3 and Ua = 0.7 cosmology is used. We use a 
Salpeter ( [Salpeter[[l955[ ) IME and all magnitudes are in the 
AB system, unless noted otherwise. 


2. SAMPLE AND SPECTROSCOPIC OBSERVATIONS 
2.1. The luminous Lya candidates at z = 6.6 



in CQSMQS/UltraVISTA (PI: M. Quchi) and - 1 deg^ in 
UDS/SXDE (PI: M. Quchi) fields in order to obtain the largest 
sample of luminous hya emitters at the epoch of re-ionisation. 

Qut of the 135 hya candidates found in Matthee et al. 
([2015j), we discover two very bright hya candidates in the 
CQSMQS/UltraVISTA field: ‘CR7’ (CQSMQS Redshift 7) 
and MASQS/^ MASQSA is particularly compact (0.7 "). 
while CR7 is extended (^ 3")- We show the location of the 
hya emitters within the CQSMQS field footprint in Pigure[^ 
in which the size of the symbols scales with luminosity. We 
also show their properties in Tabled 
Thumbnails in various wavelengths ranging from observed 
optical to observed mid-infrared are shown in Eigurej^ Both 
candidates show very high rest-frame hya EW^ in excess 
of > 200 A. By taking advantage of the wealth of data in the 
CQSMQS/UltraVISTA field (e.g.[Capak et al.|2007l[Scoville 


[et al.|2007l[Ilbert et al.|2009t [McCracken et al.|2012[ i, we ob 


^ The nickname MA SQSA consists of the initials of the first three authors 
of[Matthee et aL[12015} 

^ EWs computed by using either z or Y lead to results in excess of 200 A. 
We also present EWs computed based on Y band and from our spectroscopic 
follow-up in Table[^ 
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TABLE 1 


Sources and observation log 

Sourc^ 

R.A. 

(J2000) 

Dec. 

(J2000) 

Int. time VLT|^[Keck] 
(ks/pixel) 

Dates of observations 

Features detected in spectrunj^ 

CR7 

CR7 

CR7 

MASOSA 

10 00 58.005 
10 00 58.005 

10 00 58.005 
10 0124.801 

+ 014815.251 
+ 014815.251 
+ 014815.251 
+ 02 3145.340 

VIS: 8.1, [5.4] 

NIR: 9.9 

SINFONI: 13.5 

VIS: 6, [2.7] 

28 Dec 2014; 22 Jan, 15 Feb 2015 

22 Jan, 15 Feb 2015 

8, 11-13, 17 Mar 2015 and 4 Apr 2015 
29 Dec 2014; 12 Jan 2015 

916-1017A, Lya, Helll640A 
Helll640A 

Helll640A 

Lyo; 


^Observation log for the two luminous Lya emitter candidates, observed with both the VLT (using X-SHOOTER and SINFONI: CR7 and FORS2: MASOSA) 
and Keck (using DEIMOS for both sources) and which we have spectroscopically conhrmed. 

•^Note that for X-SHOOTER (CR7) the VIS (visible) and NIR (near-infrared) arms provide different total exposure times (due to different read-out times), and 
we thus provide them separately in different lines. Keck/DEIMOS exposure times are presented between square brackets: [exposure time]. 

^We also note the features and lines detected in the spectra: for CR7 we detect rest-frame UV just redder of the Lyman-limit, with a clear Lya forest, but such 
flux completely disappears towards redder wavelengths, implying a very blue spectra. No continuum is detected apart from that. An extremely luminous, high 
EW and narrow Lya emission line is detected for each of the sources. For CR7, where we have NIR coverage, from both X-SHOOTER and SINFONI, we detect 
high EW Heiil640A emission fro both instruments, but no other lines. 



Fig. 1.— Projected positions on sky of all Lya candidates (green circles) 
found in the COSMOS/UltraVISTA field. The grey background points rep¬ 
resent all detected sources with the NB921 filter, highlighting the masking 
applied ( due to the presence o f artefacts caused by bright stars and noisy re¬ 
gions, seejMatthee et al.|2015|. Lya candidates are plotted with a symbol size 
proportional to their Lya lurninosity. CR7 and MASOSA are highlighted in 
red: these are the most luminous sources found in the field. Their coordinates 
are given in Table 


tain multi-band photometry for both sources. The measure¬ 
ments are given in Table We re-measure the 3.6 yum and 
4.5 yum photometry for CR7, in order to remove contamina¬ 
tion from a nearby source. Such contamination is at the level 
of 10-20%, and is added in quadrature to the photometry er¬ 
rors. 

We find that CR7 has very clear detections in the near- 
infrared and mid-infrared (Figu re |^, showing a robust 
Lyman-break ( [Steidel et al.|1996 ). CR7 is detected in IRAC , 
with colors as expected for a z ^6.6 source ( |Smit et al.|2014| ), 
likely due to contribution from strong nebular lines with EWs 
in excess of a few 100 A in the rest-frame optical (see Fig¬ 
ure I3- Because of the detections in the NIR and MIR, the 
rest-frame UV counterpart of CR7 was already identified as a 
z ^ 6 —7 Lyman-break candidate ( [Bowler et k.|20 12 2014] ). 
However, because of its very uncommon NIR colors (i.e. ex¬ 
cess in J relative to V, 77 and K), the clear IRAC detec¬ 
tions, and, particularly, without the NB921 data (Figure 0, 
it was classed as an unreliable candidate, possibly a potential 


interloper or cool star. MASOSA has a clear detection in the 
narrow-band and is weakly detected in 2 ;, but the z band detec¬ 
tion can be fully explained by Lya. It is not detected at > la 
in the NIR (J > 25.7, H > 24.5, K > 24.4), although a 
very weak signal is visible in the thumbnails (Figure [^. This 
indicates that the Lya EW is very high and highlights that the 
Lyman-break selection can easily miss such sources, even if 
they are extremely bright and compact in Lya, as MASOSA 
is not detected at th e current depth of the Ultra VISTA survey 
( [Bowler et al.|2Q14| ). 

2.2. Spectroscopic observations and data reduction 

Spectroscopic observations were made with the Very 
Large Telescope (VLTfl using X-SHOOTER and SINFONI 
(for CR7) and FORS2 (for MASOSA). The choice of X- 
SHOOTER for CR7 was due to the fact that it was detected in 
the NIR and showed evidence for excess in the J band, likely 
indicating strong emission lines - SINFONI was used to con¬ 
firm the results and avoid potential biases in the slit choice 
and inclination. Both sources were observed with DEIMOS 
on the Keck II telescope (see Table[T]) as well. Spectra for both 
sources, obtained with both the VLT and Keck, are shown in 
Figure including the spectra obtained by combining both 
data-sets. 


2.2.1. DEIMOS/Keck observations 


DEIMOS/Keck observations targeted both ‘CR7’ and ‘MA¬ 
SOSA’ in two different masks and two different nights. Ob¬ 
servations were conducted on 28 and 29 December 2014. The 
seeing was ^ 0.5" on the first night, when we observed 
‘CR7’, and ^ 0.7" on the second night, when we observed 
‘MASOSA’. Observations were done under clear conditions 
with midpoint airmass of < 1.1 for both sources. We used 


a central wavelength of 7200 A and the 6001 grating, with a 
resolution of 0.65Apix“^, which allowed us to probe from 
4550 A to 9850 A. We used the 0.75" slit. 

For CR7, we obtained 4 individual exposures of 1.2 ks and 
one exposure of 0.6ks, resulting in a total of 5.4 ks. For MA¬ 
SOSA, we obtained a total of 2.7 ks. A strong, extended and 
asymmetric line is clearly seen in every individual 1.2 ks ex¬ 
posure prior to any data reduction. 

We reduced the data using the DEIM OS spec2d pipeline 
([Cooper et al. 120121 [Newman et al.[2013[). The observed spec- 
tra were flat-nelded, cosmic-ray-removed, sky-subtracted and 


^ Observations conducted under ESO DDT programs 294.A-5018 and 
294.A-5039; PI: D. Sobral. 
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Fig. 2.— Thumbnails of both luminous Lyo; emitters in the optical to MIR from left to right. Each thumbnail is 8 x S", corresponding to ~ 44 x 44kpc at 
2 ; ~ 6.6. Note that while for MASOSA the Lya emission line is detected by the NB921 filter at full transmission, for CR7 the Lya is only detected at ~ 50% 
transmission. Therefore, the NB921 only captures ~ 50% of the Lya flux: the observed flux coming from the source is ~ 2x larger. 


wavelength-calibrated on a slit-by-slit basis. We used stan¬ 
dard Kr, Xe, Ar and Ne arc lamps for wavelength solution 
and calibration. No dithering pattern was used for sky sub¬ 
traction. The pipeline also generates ID spectrum extraction 
from the reduc ed 2D per slit , and we use the optimal extrac¬ 
tion algorithm ( |Horne|p^86| ). This extraction creates a one¬ 
dimensional spectrum of the target, containing the summed 
flux at each wavelength in an optimised window. We also ex¬ 
tract the spectrum of both sources with varying apertures and 
at various positions, in order to take advantage of the fact that 
the sources are clearly spatially resolved. The flnal spectrum 
is shown in Figure 


2.2.2. FORS2/VLT observations 


FORS2/VLT ( [Appenzeller et al.||1998| ) observations tar¬ 
geted ‘MASOSA and were obtained on 12 January and 11 
February 2015. The seeing was 0.7" and observations were 
done under clear conditions. We obtained individual expo¬ 
sures of 1 ks and applied 3 different offsets along the slit. In 
total, we obtained 6 ks. We used the OG590-f 32 Alter together 
with the GRIS300I-I-11 Grism (1.62 Apix“^) with the 1" slit. 
Lya is clearly seen in each individual exposure of 1 ks. 

We use the ESO FORS2 pipeline to reduce the data, along 
with a combination of Python scripts to combine the 2D and 
extract the ID. The steps implements follow a similar proce¬ 
dure to that used for DEIMOS. 


2.2.3. X-SHOOTERNLT observations 


Our X-SHOOTER/VLT ( |Vemet et~ar |201 1| ) observations 
targeted ‘CR7’ and were obtained on 22 January 2015 and 
15 Eebruary 2015. The seeing varied between 0.8" and 0.9" 
and observations were done under clear conditions. We ob¬ 
tained individual exposures of 0.27 ks for the optical arm, 
while for NIR we used individual exposures of 0.11 ks. We 
nodded from an A to a B position, including a small jitter box 
in order to always expose on different pixels. We used 0.9 " 
slits for both the optical and near-infrared arms (resolution of 
R ^ 7500 and R ~ 5300, for the optical and near-infrared 
arms, respectively). In total, for the X-SHOOTER data, we 
obtained 8.1 ks in the optical and 9.9 ks in the NIR. The dif¬ 
ferences are driven by the slower read-out time in the optical 
CCD compared to the NIR detector. 

We use the ESO X-SHOOTER pipeline to fully reduce the 
visible (optical) and NIR spectra separately. The flnal spec¬ 
trum is shown in Eigurej^ 


We have also observed CR7 w ith the SINEONI ( [Eisenhau^ 


et al.|200^ [Bonnet et al.|2004| ) integral held unit on the VET 


on 8,11-13, 17 March 2015 and 4 April 2015. The seeing var¬ 
ied between 0.6" and 0.9" (median: 0.77") in the J band and 
observations were done under clear conditions. We used the 
non-adaptive optics mode (spaxel size: 0.25", held of view 
of 8 X 8") with the J band grism {R ^ 2000) and individual 
exposure times of 0.3 ks. We took advantage of the relatively 
large spatial coverage to conduct our observations with a jitter 
box of 2" (9 different positions for each set of 2.7 ks observa¬ 
tions). We obtained 45 exposures of 0.3 ks each, resulting in 
a total exposure time of 13.5 ks. 

We use the SINEONI pipeline (v2.5.2) in order to reduce 
the data. The SINEONI pipeline dark subtracts, extracts the 
slices, wavelength calibrates, flat-flelds and sky-subtracts the 
data. Elux calibration for each observation was carried out us¬ 
ing standard star observations which were taken immediately 
before or after the science frames. A flnal stacked data-cube 
is produced by co-adding reduced data from all the observa¬ 
tions. The collapsed data-cube does not result in any contin¬ 
uum detection, as expected given the faint J flux. We extract 
the ID spectrum using an aperture of 1". 

3. MEASUREMENTS AND SED FITTING 


3.1. Redshifts 

In both the VET (EORS2 or X-SHOOTER) and Keck 
(DEIMOS) spectra for the two targets, we detect the very 
strong Lya line (Eigurej^ in emission, and no continuum ei¬ 
ther directly red-ward or blue-ward of Lya. The very clear 
asymmetric profiles leave no doubts about them being Lya 
and about the secure redshift (Eigurej^. Particularly for CR7, 
the high S/N > 150cr (combined Ke^ and VET) at Lya, de¬ 
spite the very modest exposure time for such a high-redshift 
galaxy, clearly reveals this source is unique. 

Based on Lya, we obtain redshifts of z = 6.604 for CR'0 
and z = 6.541 for MASOSA. The redshift determination 
yields the same answer for both our data-sets: X-SHOOTER 
and DEIMOS, for CR7 and EORS2 and DEIMOS, for MA¬ 
SOSA (see Eigurej^ which shows the agreement). It is worth 
noting that for CR7 we And that the Lya emission line is de¬ 
tected in a lower transmission region of the NB921 Alter pro¬ 
file (50% of peak transmission). Therefore, the Lya luminos¬ 
ity of CR7 is higher than estimated from the NB921 photom¬ 
etry, making the source even more luminous than thought. 


3.2. Spectral line measurements 


2.2.4. SINFONI/VLT observations 


CR7 has a redshift of 2 : = 6.600 based on HeiiI640 A (see ^3.3| . 
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Fig. 3.— Left: ‘CR7’ 1-D and 2-D optical spectra, showing the strong and clear Lya emission line. We also show the NB921 filter profile which was used to 
select the source. Note that Lyct is detected at the wing of the NB921 filter. Thus, while the NB921 photometry already implied the source was very luminous, 
its true luminosity was still underestimated by a factor of two. We show both our Keck/DEIMOS and VLT/X-SHOOTER spectra, which show perfect agreement, 
but with X-SHOOTER providing an even higher spectral resolution, while the DEIMOS spectrum gives an even higher S/N. Right: ‘MASOSA’ 1-D and 2-D 
optical spectra (FORS2), showing the strong and clear Lyo: emission line. We also show the NB921 filter profile which was used to select the source. We show 
both the VLT/FORS2 and Keck/DEIMOS spectra, showing they agree very well. The DEIMOS spectrum provides higher resolution, but both clearly reveal the 
asymmetry of the line, confirming it as Lya without any doubt. 


By fitting a Gaussian profile to the emission lines, we mea¬ 
sure the EW (lower limits, as no continuum is detected) and 
FWHM. Emission line fiuxes are obtained by using NB92 1 
and Y photometry (similarly to e.g. Ouchi et al.|2009 2Q13| ), 
in combination with the NB921 filter profile and the appropri- 
ate redshift. We also check that the integrated emission line 
(without any assumption on the fitting function) provides re¬ 
sults which are fully consistent. 

Eor MASOSA, we find no other line in the optical spec¬ 
trum, and also find no continuum at any wavelength probed 
(see Figure [^. For CR7, we find no continuum either di¬ 
rectly blue-ward or red-ward of Fya in the optical spectrum 
(both in X-SHOOTER and DEIMOS; Figurej^. However, we 
make a continuum detection (spatially very compact) in the 
rest-frame 916-1017A for CR7 (rest-frame Lyman-Werner 
photons), with clear absorption features corresponding to the 
Lya forest. The reddest wavelength for which we can see 
continuum directly from the spectra corresponds to Lya at 
2 ; = 5.3. For higher redshifts, the flux is consistent with 
zero for our spectra. This clear continuum detection at wave¬ 
lengths slightly redder than the Lyman-limit, but then disap¬ 
pearing for longer wavelengths, can be explained by a com¬ 
bination of a very blue, strong continuum (intense Lyman- 
Werner radiation) and an average increase of the neutral Hy¬ 
drogen fraction along the line of sight towards higher redshift, 
similar to the Gunn-Peterson trough observ ed in quasar spec¬ 
tra (e.g. [Becker et al.|2001t Meiksin|2005| . However, due to 
the average transmission or the IGM, the fact that even just 
a fraction of the light is able to reach us is a unique finding. 
These findings will be investigated separately in greater detail 
in a future paper, including further follow-up which will allow 
an even higher S/N. 


reveals a clear J band excess (0.4 ± 0.13 mag brighter than 
expected from F, H and K\ see Table [^, which could poten¬ 
tially be explained by strong emission lines (e.g. CiV 1549 A, 
Hen 1640 A, Oiii] 1661 A, Oiii] 1666 A, Niii] 1750 A). 

We mask all regions for which the error spectrum is too 
large (> 1.5 x the error on OH line free regions), including 
the strongest OH lines. We then inspect the spectrum for any 
emission lines. We find an emission line at 12464 A (see Fig¬ 
ure 4). We find no other emission lines in the spectrum (Fig¬ 
ure 4). The emission line found, at z = 6.60C[^ corresponds 
to 1640 A, and thus we associate it with Hell. Given the line 
flux (4.1 ± 0.7 X 10“^^ erg s“^ cm“^), and the level of con¬ 
tinuum estimated from e.g. Y and H bands, the line flux we 
measure is sufficient to explain the J band excess. It also 
means that we detect Hell with a high rest-frame EW (> 20A) 
with our X-SHOOTER spectra, consistent with the neces¬ 
sary rest-frame EW in order to produce the excess in the J 
band (~ 80 A, reliably estimated from photometry). We find 
that the Hell 1640 A emission line is narrow (130 ± 30 km s“^ 
FWHM), and detected at 6 cr in our X-SHOOTER data. Re¬ 
sults are presented in Table 

In order to further confirm the reality, strength and flux of 
the Hell 1640 A line, we also observed CR7 with SINFONI 
on the VET. SINFONI reveals an emission line at the cen¬ 
tral position of CR7 (at the peak of Lya emission), spatially 
compact (unresolved) and found at 12464 A, thus matching 
our X-SHOOTER results. We co-add the X-SHOOTER and 
SINFONI spectra and show the results in Figure When 
co-adding the spectra, we normalise both spectra at the peak 
value of Hell 1640 A and explicitly mask the strong OH lines. 


3.3. NIR Spectra ofCR7: Hell and no other lines 

We explore our X-SHOOTER NIR spectra to look for any 
other emission lines in the spectrum of CR7. The photometry 


The redshift measured from the Hell 1640 A emission line implies a 
small positive velocity offset between the peak of Lyo; and Hell 1640 A of 
+160kms“^, i.e., the Lya peak is redshifted by +160kms“^ in respect to 
Hen. 
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TABLE 2 

A SUMMARY OF OUR RESULTS FOR CR7 AND MASOSA AND 
COMPARISON TO HiMIKO. 


Measuremenj^ 


CR7 

MASOSA 

Himiko 


Zspec Lyo; 


D.0U4_Q QQ3 

6.541 ± 0.001 

6.54 


/3 UV slope 


-2.3 ±0.08 

— 

-2.0 ±0.57 


Lya (EWHM, km s 

-1) 

266 ± 15 

386 ± 30 

251 ±21 


Lya (EWo,o6«,ltJ A) 

211 ±20 

> 206 

78±8 


Lya (EWq gpgc? 

> 230 

> 200 

— 


Lya (LogioL, ergs 


43.93 ± 0.05 

43.38 ± 0.06 

43.40 ± O.Of 

Lya/NV 


> 70 

— 

— 


Heii/Lya 


0.23 ±0.10 

— 

— 


Hen (EWo, A) 


80 ± 20 0 20) 

— 

— 


Hen (EWHM, km s 


130 ± 30 

— 

— 


Heii/OIII]1663 


> 3 

— 

— 


Heii/CIII]1908 


> 2.5 

— 

— 


Photometry 


CR7 

MASOSA 

Himiko 


z 


25.35 ±0.20 

26.28 ±0.37 

25.86 ± 0 . 2 c 

i 

NB921 2" aperture 

23.70 ±0.04 

23.84 ±0.04 

23.95 ± O.Oi 

? 

NB921 MAG-AUTO 

23.24 ±0.03 

23.81 ±0.03 

23.55 ± 0.05’ 

r 

Y 


24.92 ±0.13 

> 26.35 

25.0 ±0.3^ 

1 

J 


24.62 ±0.10 

> 26.15 

25.03 ±0.2^ 

1 

H 


25.08 ±0.14 

> 25.85 

25.5 ± 0.35*' 


K 


25.15 ±0.15 

> 25.65 

24.77 ±0.29' 

i 

3.6/j,m 


23.86 ±0.17 

> 25.6 

23.69 ± 0.09' 

i 

4.5/im 


24.52 ±0.61 

> 25.1 

24.28 ±0.19' 

i 


^These include both the spectroscopic measurements, but also photometry. 
In order to provide an easy comparison, we also provide the measu rements 
for Himiko, the o ther luminous sour ce in thejMatthee et al.|(|2015) sample, 
fully presented in|Ouchi et al.U2013). Note that the intrinsic Lyo; properties 
may be significantly ditterent from the observed ones, if there is significant 
absorption by the galaxies themselves. 

Using Y band photometry to estimate the continuum. 

^Re-computed by using Y band to estimate the continuum, in order to match 
our calculation for C R7 and MASOSA. 

^Measurement fr om|Ouchi et al. j2013) 

^Measurement bv| Maiihee ei al. (2 015p b ased on own UDS reduction. 

^Measurement from lBowler et aiTj2U14) since the Y and H data has been 
significantly improved with the availability of new VISTA and UKIDSS data. 

Our Hell 1640 A emission line implies a small velocity off¬ 
set between the peak of Lyct and Hell 1640 A of - 1 -160 km s“^ 
(i.e., the Lya peak is redshifted by -fl60kms“^ in respect to 
Hell, which could also be interpreted as an outflow). This 
means that, not surprisingly, we are only detecting the red 
wing of the hya line, while both the blue wing and any po¬ 
tential Lya component with < +160kms“^ is being likely 
absorbed and cannot be observed. This may imply that the in¬ 
trinsic Lya flux or luminosity and the EW will be even larger 
than measured. 


4. DISCOVERY OF THE MOST LUMINOUS LYa EMITTERS 
4.1. MASOSA 

MASOSA is particularly compact (0.7 " in diameter, cor¬ 
responding to 3.8 kpc diameter at 2 : = 6.541). It is similar to 
sources now found by MUSE (e.g. [Bacon et al. Karman 

jet al.|2015] ), that are completely undetected in the continuum 
and thus similar to typical Lya emitters that have been found 
over the last years with e.g. Subaru/Suprime-cam. However, 
MASOSA is extremely bright in Lya and thus separates it 
from the very faint sources found with MUSE in the HDE 
South and from other Lya emitters. MASOSA is undetected 
in all continuum bands at all wavelengths and even the weak 
detection in z' can be fully explained by the luminous Lya 
line and little to no continuum. The estimated rest-frame EW 


of the Lya line from spectroscopy is very high (> 200 A), thus 
implying a likely v ery metal-poor stellar population ( jMalho- 
tra & Rhoads [2002 j). Its nature is likely similar to other metal- 


s)2002| ). 

litters (|? 


poor Lya emitters ( jNagao et al.|2008[|0no et al.|2010| 2012] ). 

-Given the very high equivalent width (EW) and no continuum 
-detection, mAsOSA is likely extremely young, metal-poor 
and likely contains low stellar mass (< 10^ Mq). 

The Lya emission line shows a EWHM of 386 ± 30 km s“^, 
and with tentative evidence for t wo components in the Lya 
line (similar to what is found by Quchi et al. 2010| ), p oten¬ 
tially indicating a m erger or radiative transfer etfects (e.g. jVer-j 
jhamme et al.||2006| ). However, with the current ground-base 
imaging, and the faintness in the continuum, it is not possible 
to conclude anything about the potential merging nature of 
the source - only HST follow-up can investigate this. How¬ 
ever, what is already clear is that MASOSA provides a new 
-example of a relativel y compact Lya emitter w ith a similar 
-luminosity to Himiko ( Quchi et'n^|2009||2Q13 ), but with a 
much more extreme EW, revealing that such high luminos¬ 
ity sources present diverse properties and may have a diverse 
nature. 


4.2. CR7 

CR7 clearly stands out as the most luminous Lya emitter 
at ^ ^ 7, with a luminosity oiL^ya = ergs“^, 

3 X more luminous than any known Lya emitter within the 
epoch of re-ionisation. It also presents a very high rest-frame 
EW of > 200 A, and thus a likely intrinsic EW which is even 
higher because of absorption by the ISM. Our measurements 
are presented in Table 

CR7 is spatially extended: 3" in diameter, corresponding 
to ^ 15 kpc at 2 ; = 6.6, as seen from its narrow-band image, 
but also in the spectra. Both the X-SHOOTER and DEIMOS 
data confirm its spatial extent, with both of them agreeing per¬ 
fectly on the redshift, extent and EWHM. While the Lya line 
profile is narrow (EWHM« 270kms“^), and particularly for 
such high luminosity, we see evidence for potentially 2 or 3 
components (double peaked Lya emission and a redshifted 
component towards the South of the source) and/or signs of 
absorption (see Eigurej^, which indicate a complex dynami¬ 
cal structure. It may also mean that the actual intrinsic EWHM 
is even narrower. However, such tentative evidence requires 
confirmation with further spectroscopy obtained over differ¬ 
ent angles, and particularly by exploring deep imaging with 
high enough spatial resolution with e.g. HST. 


4.3. Comparison with Himiko 

We find that CR7 may be seen as similar to Himiko (but 
much brighter in Lya and much higher EW) due to both 
sources presenting a spatial extent of about 3" in diameter. 
Both could therefore be tentatively classed as Lya “blobs” 
(e.g. Steidel et~ST]|2000t [Matsuda et al.||2004} jSteidel et~ST 


20111. However, we note that Himiko is detected at peak 


transmission in the NB and t he NB imaging in w hich it is 
detected is 1 mag deeper (see [Matthee et al.||2015| ) than the 
imaging used for the discovery of CR7. While the CR7 Lya 
line profile is very narrow, it consists of 2 or potentially 3 
components, which may indicate that the source is a double or 
triple mer ger, likely similar to Himiko in that respect as well 
(see ^5.4 which shows this is very likely the case for CR7). 
However, a simpler explanation is radiation transfer , which 
can easily cause such bumps (e.g. jVanzella et al.|2010| ). There 
are other similarities to Himiko, including: detections in NIR 
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Fig. 4.— Left: Our X-SHOOTER NIR spectrum of CR7, revealing a significant detection of the Heiil640 A emission line. We show both the data at full 
resolution and binning in wavelength (with a resolution of 0.4 A). We also show the combined X-SHOOTER and SINFONI data and also show the sky spectrum. 
We note that sky lines were explicitly masked, but some residuals are still visible, including those of two OH lines just redder of the Hell emission line, which 
show up as a slight flux increase and another one bluer of Hell, which has been slightly over-subtracted. We note that no continuum is found in the NIR 
spectra, and that we would require significantly deep er observations in order to detect it simply based on the NIR photometry. Right: We investigate the NIR 
X-SHOOTER spectrum for other emission lines (e.g.|Stark et al.|2014|. We do not find any other emission line apart from Hell, but we show one of the lines 
that should be stronger in our spectrum (typically ~ 2 x stronger than Hell, thus Heii/Ciii]~ 0.5). We use our Hell line detection and show it at the position of 
the Cm] doublet for a typical line ratio of Heii/Ciii] of 0.5, but also a line ratio of ~ 1 and ~ 2. Our data allow us to place a limit of Heii/Ciii]> 2.5. We also 
investigate the presence of NV1240, NIV1487, CIV1549, Oiii]1661, Oiii]1666, Niii]1750: all these are undetected. 


and a blue IRAC color. CR7 is however a factor ^ 3 brighter 
in Lyct emission and has an excess in J-band (attributed to 
Hell emission). CR7 is also bluer = —2.3 ± 0.08, ei- 
ther using Y — H or H — K, following equation 1 of Ono 
et al.|201Q| in the rest-frame UV when compared to HimlEo 
(which shows P ^ —2.0, but note that Himiko shows a red 
colour from H io K which would imply {3 ^ 0.2 if those 
bands are used). Furthermore, while for CR7 we see some 
rest-frame UV light just redder of the Lyman-limit, corre¬ 
sponding to Lyman-Werner radiation (very compact, coinci¬ 
dent with the peak of Lya emissio n and the HST det ection), 
this is not seen at all for Himiko ( |Zabl et J^|2015| ). Also, 
while CR7 shows a strong HeHl6 40 A emission line , no Hell 
is detected in Himiko, even though [Zabl et al.| ( |2015| ) obtained 
very deep X-SHOOTER data. 

MASOSA is quite different. While it has the highest Lya 
peak brightness, it is not extended and not detected in NIR 
or IRAC. Therefore, MASOSA provides also a new class of 
sources at the epoch of re-ionisation: as luminous as Himiko, 
but very compact and with no significant rest-frame UV or 
rest-frame optical detection at the current UltraVISTA depth. 

5. SED EITTING AND MODEL ASSUMPTIONS 

To interpret the ph otometry/SED of CR7 w e explo it the 
SED-fitting code of |Schaerer & de Barros| ( |2009| ) and 
[Schaerer & de Barros| ( |2010| ), which is based on a ver¬ 
sion of iho Hyperz photometric redshift code of [Bolzonella 
jet al.| ( |2Q00| ), modified to take nebular emission into account. 
We have explored a variety of spectral template s including 
those from the GALAXEV synthesis models ofjBruzual & 
jCharlotj (|2003|), covering different metallicities (solar, Zq, 
to 1/200 Zq) and star formation histories (bursts, exponen¬ 
tially declining, exponentially rising). A standard IMF with 
a Salpeter slope from 0.1 to 100 Mq is assumed. We refer 
to these models as “standard”/”enriched” SED fits or “stan- 
dard”/”enriched” models throughout this paper. 


In addition, we also use synthetic spectra from metal- 
free (PopIII) stellar populations a ssuming different IMF s 
(Salpeter, top-heavy), taken from jSchaererj ( [2002 2003| ). 
Constant SFR or bursts are explored in this case. We also 
explore SEDs from “composite” stellar populations, showing 
a superposition/mix of PopIII and more normal populations. 

Nebular emission from continuum processes and emis- 
sion lines are added to the sp ectra predicted as described in 
[Schaerer & de Barrosj ( [2009| ). Nebular emission from con¬ 
tinuum processes and emission lines are proportional to the 
Lyman continuum photon production. Whereas many emis¬ 
sion lines are included in ge neral, we only i nclude H and He 
lines for the PopIII case (cf. [Schaerer|2003[ ). The int ergalac 


tic me dium (IGM) is treated with the prescription of [Madau 
( 1995[). Attenuation by dust is described by the Calzkti law 
( Calzetti et al.|2000] ). 


5.1. CR7: SED fitting with a normal population 

Part of the photometry of CR7 is explained relatively well 
with “standard” models, as illustrated in Figure[^(left panel). 
The shortcomings of these fits is that they cannot account for 
the relative excess in the J band with respect to V, 77, and 77, 
and that the Eya emission is not strong enough to reproduce 
the entire flux observed in the NB921 filter. Both of these 
shortcomings thus relate to the presence of the strong emis¬ 
sion lines (Lya and Hell) observed in spectroscopy and also 
affecting the broad-band photometry. 

The typical physical parameters derived from these SED 
fits which only include “normal” stellar populations indicate 
a stellar mass ^ 2 x 10^^ Mq, SFR ^ 25 Mq yr“^, and a 
fairly old age (~ 700 Myr; the Universe is ^ 800 Myr old at 
z = 6.6). The SED fits and the derived parameters do not vary 
much for different star-formation histories (SFHs): both for 
exponentially declining and rising cases, the fits prefer long 
timescales (i.e. slowly varying SFHs). Whereas for exponen- 
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Fig. 5.— Left: The SED of ‘CR7’, from observed optical (rest-frame FUV) to observed MIR (rest-frame optical) and the best fit with a with a normal stellar 
population (not including PopIII stars). Red crosses indicate the flux predicted for each broad-band hlter for the best ht. The ht fails to reproduce the strong Lya 
emission line and also the excess in J band, due to Hell emission. Moreover, and even though the fit is unable to reproduce all the information available for the 
source, it requires an age of 700 Myr (the age of the Universe is 800 Myr at 2: = 6.6). In this case, the galaxy would have a SFR of ~ 25 Mq yr“^ and a stellar 
mass of ~ 10^°-^ Mq. This, of course, is not able to explain the strong Lya and the strong Hell emission line. Right: Same observed SED of CR7 as in the 
left panel plus HST photometry for clumps B-i-C (magenta triangles). The black line shows a ht with a pure PopIII SED to the rest-frame UV part; the magenta 
line the SED of an old simple stellar population with 1/5 solar metallicity which matches the hux from clumps B-i-C; the green line shows the predicted SED 
summing the two populations after rescaling the PopIII SED by a factor 0.8. The composite SED reproduces well the observed photometry. Although there is 
a tension between the strength of the Heii line and nebular continuum emission (cf. text), a PopIII contribution is required to explain the Hell A1640 line and 
the corresponding excess in the J-band. Although Hell is very strong, we hnd no evidence for any other emission lines that would be characteristic of an AGN. 
Furthermore, the clear IRAC detections and colors, and particularly when taken as a whole, can be fully explained by a PopIII population, while normal stellar 
population or AGN is simply not able to. 


tially declining SFHs and for constant SFR the best-fit atten 
nation is negligible {Ay = 0), a higher attenuation is needed 
for rising SFHs, as expected (cf. Schaerer & Pell6|2005[ Fin- 
lator et al.|2007| ). Depending on the assumed metallicity, this 
may reach from Ay = 0.5 (for 1/5 solar) to 0.25 (for 1/200 
Zq). The corresponding SFR is ^ 30 — 4OM0yr“^, and 
the predicted IR luminosity ranges from log(I/iR/ Lq) = 10.4 
to 11.2, assuming energy conservation. The typical specific 
SFR, sSFR, obtained from these fits is sSFR ^ 1.2 — 1.9 
Gyr“^, as expected for a “mature” stellar pop ulation with 
SFH close to constant (c.f. [Gonzalez et al.|2014| . 

5.2. CR7: SED fitting with contribution from PopIII stars 

The presence of strong Eya and Hell emission lines, plus 
the absence of other UV metal emission lines (cf. above), may 
be due to exceptionally hot stars with a strong and har d ion¬ 
izing fiux, resembling that expect ed for PopIII stars (cf. Tum-j 
linson et al. [200 1 [ |Schaerer|2002| ) . A fit with Pc^III templates 
(from|Schaerer|20G 

black line). PopIII 
at 3.6 yum (for 2 ; = 6.6) there are strong Hell and Hel lines, 
apart from H/5, and that these Helium emission lines should 
have fluxes com parable to that of H/3 in the case of PopIII 
( |Schaerer|2002| ). Specifically, 3.6 pm should be contaminated 
by HeH4686, Hei4471 and Hei5016. For 4.5 pm, apart from 
Hof, the Hei5876 should be detected and could be compara¬ 
ble to H/3. The general features of our fits with PopIII tem¬ 
plates (considering bursts or constant SFR, as well as different 
IMFs) are the following: 

• The UV rest-frame part of the SED is very well re¬ 


et ai.izuuip^cnaererizuuzp. A nt witn uopm templates 
Sctiaerer|2002||2003 I is shown in Figure [5| (right panel: 
line). PopIII models ( |Schaerer||2002 2(^3] ) show that 
re are strong Hell and 


produced, allowing naturally for a “boost” of the J- 
band fiux due to the presence of the Hell line and for 
a stronger Lya line due to a higher Lyman continuum 
fiux (for the same UV fiux). The exact strength of these 
emission lines is very sensitive to the “details” of the 
population, such as the upper end of the IMF, the age or 
star formati on history (see more details in e.g. [Raiter 
et al.j 2010| ). The fits with PopIII also reproduce all 
other NIR detections: Y, H and K (see Figure [^. 

• Populations with very strong line-emission also have 
a strong nebular continuum emission red-ward of hya 
and increasing towards the Baimer limit (see e.g. 
|Schaerer||20Q2j ). The observations (H and K band 
photometry) do not permit a much stronger contribu¬ 
tion from the nebular continuum, as that would mean 
an increasing redder H-K colour, which is not ob¬ 
served. This limits the maximum strength of the pre¬ 
dicted emission lines, except if the two emission pro¬ 
cesses (recombination line emission and nebular con¬ 
tinuum emission, which is due to two-photon and free- 
bound emission) could be decoupled, and emission 
lines could be increased without significantly increas¬ 
ing the nebular continuum emission. 

• All the available PopIII templates fitting the (rest-)UV 
part of the spectrum, predict a relatively low fiux in the 
IRAC bands with or without accounting for emission 
lines. Therefore, a pure metal-free population does not 
seem to be able to completely reproduce the observed 
rest frame UV-optical SED of this source. In any case. 
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Fig. 6.— Left: CR7 with the NB921 filter/Suprime-cam imaging on Subaru, showing the extent of the Lya but note that NB921 detects Lya at only 50% 
transmission. Middle: HST imaging in Y J, revealing that CR7 clearly splits into 3 different components which we name A, B and C. Right: HST imaging in 
H, again revealing the 3 different components in CR7. We find that component A fully dominates the rest-frame UV and is coincident with the peak of Lya 
emission and the location at which we detect strong Heiil640A emission. Clumps B and C are much redder, and fully consistent with significantly contributing 
to the IRAC photometry. Note that because of the colours of the B and C clumps, they completely dominate the mass of the system, and thus the actual mass 
centre of the system would be located between C and B, and significantly away from A. This is fully consistent with a scenario in with PopIII star formation is 
propagated in a wave from the central position towards the outskirts. 


a PopIII only explanation would not seem very likely. 
Therefore, a metal-free population alone (without de¬ 
coupling between recombination line emission and neb¬ 
ular continuum emission) is not able to reproduce the 
observed rest-frame UV-optical SED of this source. 

5.3. PopIII and a more chemically evolved stellar population 

As a consequence of our findings, we are led to consider 
a hybrid SED consisting of two populations (the source may 
well be a merger, and these components may well be com¬ 
pletely separated, avoiding pollution by metals into the po¬ 
tentially metal-free region): a young metal-free stellar popu¬ 
lation and a more chemically evolved population. This can be 
fully confirmed by using HST and appropriate filters that can 
easily isolate Lya, the rest-frame UV, and Hell. 

A superposition of a young PopIII component dominating 
in the UV and an older population of “normal” metallicity (in 
this case, 0.2 Zq) dominating the rest-frame optical flux of 
CR7 is shown in the right panel of Eigure In practice, we 
add 80% of a metal-free simple stellar population with an age 
of 16 Myr (shown by the black line; although younger ages of 
< 5 Myr are preferred to produce the strongest Hell emission) 
to a 360 Myr old burst of 1/5 Zq (magenta curve), giving the 
total fiux shown in green. As can be seen, an “old” population 
of - 1.6 X 10^^ Mq can make up for the missing rest-frame 
optical fiux, whereas a young PopIII burst can dominate the 
UV and the emission lines. There is some tension/uncertainty 
in the age or age spread of the metal-free population, as very 
young ages (< 5 Myr) are preferred to produce the strongest 
Hell emission, whereas slightly older ages are preferred to 
avoid too strong nebular continuum emission (cf. above). Eor 
indication, the mass of the metal-free component would be 
1.4 X 10^ Mq for a Salpeter IMP from 1 to 500 Mq, i.e. — 9 
% of mass of the old population. However, significantly less 
mass could be needed if the PopIII IMP was fiat or top-heavy, 
lacking e.g. completely low mass stars (stars below 10 M0). 
This could mean that — 10^ Mq of PopIII stars would be 
needed in order to fully explain the fiux for an IMP peaking at 
— 60 Mq, or even less if the IMP peaks at even higher masses. 
This reveals that the presence of a young, metal-free popula¬ 
tion, forming for example in a yet un-polluted region of the 
galaxy, in an slightly evolved galaxy at z = 6.6, could repro¬ 


duce the observed fe atures of CR7 (consistent with theoretical 
predictions from e.g.|Scannapieco et al.|2003l|Tomatore et al. 
|2007l l. 

5.4. HST imaging ofCR7 

In order for our best interpretation to be valid, CR7 would 
require to be clearly separated/resolved, with HST resolu¬ 
tion, into at least two different spatial components: one being 
dominated by a PopIII-like stellar population (dominating the 
UV light but with only a very small fraction of the mass), and 
another, redder in the UV (and fully dominating the mass), 
with fluxes similar to those shown in Pigure While this 
requires more detailed follow-up, we find that CR7 has been 
fortuitously observed and is in the PoV of previous WPC3 
observations in PllOW (broad YJ filter that also contains 
Lya) and P160W (H) of a different project (ID: 12578, PI: 
Porster-Schreiber). We explore such data in order to investi¬ 
gate the rest-frame UV morphology of CR7 and to conduct 
a first study of the rest-frame UV colours. Observations in 
PllOW and P160W were obtained for 2.6ks each. 

We show the data in Pigure including a comparison with 
out NB921 imaging. Pigure [Tjpresents a false-colour image 
combining data obtained witn NB921, PllOW and P160W. 
We find that CR7 is, beyond any doubt, split in different com¬ 
ponents (see Pigures|^and[7]), in line with our best interpre¬ 
tation of a PopIII-like stellar population which dominates the 
UV, and a redder stellar population, found to be physically 
separated by at least 5 kpc (projected). In fact, we actually 
find three different components, which we label as A, B and 
C (see Pigures[^and[^. We obtain photometry for each of the 
clumps separately, in order to quantitatively test if they could 
ex plain the UV photometry predicted for the two components 
in ^5.3 We use 0.4" apertures for components B and C and 
1" for component A, more spatially extended. 

We find that the sum of the two redder, fainter clumps 
(B+C) matches our evolved stellar population remarkably 
well (see Pigures[^and[^. Note that the photometry of the 
two redder clumps was not used to derive such fit. We find 
that the central clump (A) is the one that dominates the rest- 
frame UV light (Pigure [^. Pigure^ also shows the rest-frame 
colours of components A, B and CTWe also find that the peak 
of Hell emission line (from SINPONI, although with low spa- 
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Fig. 7.— A false colour composite of CR7 by using NB921/Suprime-cam 
imaging (Lya) and two HST/WFC3 filters: FllOW (YJ) and F160W (H). 
This shows that while component A is the one that dominates the Lya emis¬ 
sion and the rest-frame UV light, the (likely) scattered Lya emission seems to 
extend all the way to B and part of C, likely indicating a significant amount 
of gas in the system. Note that the reddest (in rest-frame UV) clump is C, 
with B having a more intermediate colour and with A being very blue in the 
rest-frame UV. 


tial S/N) and the peak of Lya (high S/N but extended) likely 
originate (within the errors) from the centre of the bright¬ 
est UV clump (A). The clumps are physically separated by 
^ 5kpc. 

The HST imaging reveals that CR7 may be either a triple 
merger (similar to Himiko), and/or a system where we are wit¬ 
nessing a PopIII star formation wave, which may have moved 
from the reddest clump (C) to the other (B) and we are ob¬ 
serving the brightest UV clump at the right time (A). We have 
direct evidence of the intense Lyman-Werner radiation (rest- 
frame 912-1000 A) from the brightest UV clump. It is there¬ 
fore possible that the other clumps have emitted as much or 
even more of such radiation a few ^ lOOMyrs before, pre¬ 
venting what is now the site of young massive stars (A) to 
form before and potentially allowing for that pocket of metal 
free gas to remain metal free. There are of course, other po¬ 
tential interpretations of our observations. In ^ we discuss 
the different potential scenarios in detail. 


6. DISCUSSION 
6.1. The nature of CR7 

CR7, with a luminosity of ergs“^ is 

^ 3 X more luminous than any known Lya em itter within the 
epoch of re-ionisation (e.g. |Ouchi et al.[20T^ . 

Our optical spectrum shows that the source is very blue to¬ 
wards the extreme ultra-violet up to the Lyman limit at 912 
A rest-frame, as we detect some faint continuum (spatially 
very compact) at rest-frame ^916 — 1017A (Lyman-Werner 
radiation). X-SHOOTER data also provides a near-infrared 
spectrum, allowing to investigate the signi ficant excess seen 
in the J band photometry from Ultra VISTA ( [McCracken et al. 
|2Q12t [Bowler et aH2014j ), indicative of emission line(s). No 
continuum is detected in the NIR spectrum. However, and de¬ 
spite the relatively low integration time, a strong Hell 1640 A 
line was found (^ 6cr), capable of explaining the excess in the 
J band (see Figurej^. Hell can only be produced if the intrin¬ 
sic extreme UV spectrum is very hard, i.e., emits a large num- 



Fig. 8.— HST imaging in Y J and H allows us the physically separate CR7 
in two very different stellar populations and show remarkable agreement with 
our best-fit composite SED derived in §5.3| While clump A (see e.g. Figure 
[tJ is very blue and dominates the rest-frame UV flux, B-i-C are red and likely 
^minate the rest-frame optical and the mass. Note that the we si mply show 
the HST data together with our best fit composite model derived in §5. 3 [ which 
was solely based on the full photometry and did not make use of any resolved 
FfXrdata. 


ber of ionising photons with energies above 54.4 eV, capable 
of ionising He completely. From our X-SHOOTER spectra, 
we place a lower limit in the rest-frame EW of the Hell line 
of > 20 A, but estimate from photometry that the line has 
EWq = 80 ± 20 A, consistent with our spectroscopic lower 
limit. The line we detect is also narrower than Lyct, with 
FWHM of 130 ± 30kms“^, as Heiil640 A does not scatter 
easily as Lya, as the line is not self-resonant. 

While in principle there are a variety of processes that 


could produce both high EW Lya and Helll640A, some of 
them are very unlikely to produce them at the luminosities 
we are observing, such as X-ray binaries or shocks. How¬ 
ever, in principle, cooling radiation could produce strong 
Lya e mission with luminosities simi lar to those measured for 
CR7. Faucher-Giguere et al.[ ( 201()| ) provides predictions of 
the total Lya cooling luminosity as a function of halo mass 
and redshift. Under the most optimistic/extreme assump¬ 
tions, it would be possible to produce a Lya luminosity of ^ 
10^^ erg s“ ^ for a dark matter halo mass of M > 5 x 10^^ Mq . 
Since such dark matter haloes should have a co-moving num¬ 
ber density of about ^10“^ Mpc“^ at 2 ; ^ 6.6, their number 
densities could potentially match the luminous Lya emitters 
that we have found. However, in the case of cooling radi¬ 
ation, the Hell emission line should be significantly weaker 
tha n what we measur e (with intrinsic Hell/Lya of 0.1 at most; 
e.g. [Yang et al.|20Q6[ ), and cooling radiation in a massive dark 
matter halo should also result in a broader Lya line than what 
we observe. 

There are, nonetheless, four main sources known to emit 
an ionising spectrum that can produce high luminosity, high 
EW nebular Lya and H ell as seen in our spectra (see also 
similar discussion in e. g. [Prescott et al.|20Q9[ Cai et al.|2011t 
[Kashikawa et al.|20T^ : 


1) strong AGN (many examples have been found, partic- 
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ularly on e.g. radio galaxies: |De Breuck et al.pOOO ), with 
typical FWHM of lines being ~ 1000 km s“^; 

2) Wolf-Rayet (WR) stars (man y case s known i n e.g . 
SDSS: Shirazi & Brinchmann|2012} or see |Erb et al.| ( |2010| ) 
for a higher redshift example), with typical FWHM of lines 
being ^ 3000 km s“^; 

3) Direct collapse black hole (DCB H), which have been 
predic ted and studied theoretically (e.g. [Agarwal et al.| 


|2015| ), although none has b een identified yet, bu t they shouL 
produce strong Hell 1640A ( Johnson et al.|2011 1; 


4) PopIII stars (e.g. |Schaerer||2003t [Raiter et~ar]|2010| ), 
which should produce high EW, narrow Hell emission lines 
(FWHM of a few - 100 km s-^). 

Many potential candidates for PopIII have been identified 
based on their c olours and/or high E W Lya, but either no Hell 
was found (e.g. Nagao et al.|2008|). Hell was found but with 
clear signatures of AGN activity ( De Breuck et al.|2000 Mat- 
suoka et ar]|2009 ), or Hell had very low EW ( [Cassata et al. 


2013 1 ). Thus, so far, not a single source has been found with 
high EW hya. Hell detection with high EW, no AGN sig¬ 
natures, no WR star signatures (e.g. P-Cygni profiles, broad 


lines, a nd many other metal line s e.g. |Brinchmann et al.|2008 
see also |Grafener & Vink|2015| ) and with no other metal lines. 

6.2. The nature ofCR7: AGN or WR stars? 

In order to test the possibility of CR7 being an AGN, we 
start by checking X-ray data. We fin d no detection in t he X- 
rays, with a limit of < lO^^ergs”^ ( Elvis et ^|2009 ). We 


also find no radio emission, although the limit is much less 
stringent than the X-ray emission. The X-SHOOTER and 
DEIMOS spectra were carefully investigated for any metal 
lines, particularly NV, Oxygen and Carbon lines (see Table 
and also Figure ffl. No such lines were found, and thus we 
place Icr upper limits on their fiuxes, to constrain the nature 
of the source, finding, e.g. Lya/NV > 70. Our Lya and Hell 
emission lines are narrow (both FWHM^ 100 — 300 km s^), 
thus excluding broad-line AGN. Narrow-line AGNs with Hell 
emission typically h ave CIII]1909/Heii^ 1.5 ± 0.5 (e.g. [Pe 
Breuck et al. |20()0| ); such line ratio would result in a strong 
C111J1909 detection in our X-SHOOTER spectrum (see Fig¬ 
ure ffl. We do not detect Cm] 1909 and obtain a strong up¬ 
per hmit of Clll]1909/Hell< 0.4 (la; see Figure]^, which 
greatly disfavours the AGN hypothesis. The strong limit on 
Lya/NV >70 also disfavours the AGN hypothesis and points 
towards very low metallicities. 

There are no indications of WR stars, due to the very narrow 
Hell line (^ lOOkrns”^, com pared to typical FWHM of ^ 
3000kms“^ for WR stars, c.f. Brinchmann et al. 2008| and 
no other metal lines. 

We note, nonetheless, that while CR7 is strongly dis¬ 
favoured as an AGN, it shows characteristics of what has been 


predicted for a direct collapse bla ck hole (e.g. [Johnson et al. 
20111 [Agarwal et al.|[MT3| [2015[). This is because it shows 


(3 = —2.3 (as predicted), no metal lines, and high luminos 
ity. The detection of the other nearby sources about ^ 5 kpc 
away (see Figure!^ are also a key prediction from Agarwal et 
al. 2015, while the relatively high observed Heii/Lya would 
also match pred ictions for a direct collapse black hole ( [John-| 
son et al.|201 1| ) However, CR7 does not show any broad line, 
as predicted by [Agarwal et al.[ ( 20131) , and the observed Lya 
and Hell luminosities are higher by about ^ 2 orders of mag- 
nitude when compared to predictions by e.g. [Johnson et al. 


([20111 for the case of ^ 1 — 5 x 10^ Mq black holes. An- 
other key distinction between a direct collapse black hole and 
stellar population(s) is X-ray emission: if it is a black hole, it 
must be emitting much more X-ray fiux than a PopIII stellar 
population, and thus would likely be detectable with Chandra 
given the high line luminosities measured. Deeper Chandra 
observations could in principle test this. 

6.3. On the hardness of the ionizing source ofCR7 

The theory of recombination lines relates, to first order, the 
ratio of He ll to Hydrogen recombination lines to the ratio 
between the ionizing photon fiux above 54 eV, (5(He+), and 
that above 13.6 eV, Q{H), the energies needed to ionize He+ 
and H respectively. For the relative intensity He ll A1640/Lya 
one thus has: 


/(1640)//(L^a) ^ 0.55 x 


Q(He+) 
Q{H) '' 


( 1 ) 


where the numerical factor depends somewhat on the elec 
tron t emperature, here taken to be = 30 kK ( [Schaerer 
2002]). The observed line ratio He ii A1640/Lya ^ 0.22 
therefore translates to Q(Re'^)/Q{H) « 0.42, which indi¬ 
cates a very hard ionizing spectrum. For metal-free stellar 
atmospheres such a hardness is only achieved in stars with 
very high effective temperatures, typically Teff > 140 kK 
(or > 100 kK for 1 a lower limit), hotter than the (already 
hot) zero-age main sequence predicted for PopIII stars which 
asymptotes to T eg 100 kK for the most massive stars (cf. 
Schaerer|2002[ ). For integrated stellar populations consisting 
of a ensemble of stars of different masses, a maximum hard¬ 
ness Q(He^ )/Q(i7) ^ 0.1 is exp ected for zero or very low 


metallicities (|Schaerer|2002 200 3|), higher than inferred from 
the observed He ll A1640/Lya ratio of CR7. This could indi¬ 


cate that only a fraction of the intrinsic Lya emission is ob¬ 
served, or that a source with a spectrum other than predicted 
by the above PopIII models (e.g. an AGN) is responsible for 
the ionization or contributing at high energies (> 54 eV). In 
fact the observed He ii A1640 equivalent width of 80 ± 20 
A is in good agreement with the maximum e quivalent width 
predicted for PopIII models ( Schaerer[[2002| ). This supports 
the explanation that ^75% of the intrinsic Lya emission may 
have escaped our observation, e.g. due to scattering by the 
IGM, internal absorption by dust, or due to a low surface 
brightness halo, proces ses which are known to affect in gen¬ 
eral Lyg emission (e.g. [Atek et al.|2008[ [Dijkstra et al.|201l[ 
Steidel et al.[|2011 ). If the intrinsic Lya emission is > 3 — 4 
times higher than observed, the hardness ratio is compatible 
with “standard” PopIII models. Future observations and de¬ 
tailed photoionization models may yield further insight on the 
properties of the ionizing source of CR7. 

6.4. CR7: A PopIII-like stellar population? 

As the AGN and WR stars hypothesis are strongly dis¬ 
favoured (although we note that a direct collapse black hole 
could still explain most of our observations), could CR7 be 
dominated by a PopIII-like stellar population? For this to 
be the case, such stellar population would have to explain 
the detection, FWHMs, EWs and fiuxes of Lya and Hell 
(including the strong excess in J band), the UV continuum 
and continuum slope {j3 = — 2.3 ± 0.08) and the IRAC de¬ 
tections which impl y ve ry hi gh E W rest-frame optical lines. 
We hav e shown in 03 and ^ 5.4 that a compos ite of PopIII 
models ( [Schaerer|2002[ 2003 [[Raiter et al.|2010| ) with a more 
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evolved stellar population can match all our observations, in¬ 
cluding spatially resolved HST data. We note that the intrinsic 
Hell/Lyg line ratio pre dicted for PopIII would be ~ 0.05 — 0.1 
(|Schaerer|2002| |2003|), but that can easily result in an observ- 
able ratio of ^ 0.2 — 0.3 if a significant fraction of the Lya 
line is absorbed/attenuated by neutral Hydrogen. Since we 
find a velocity offset between Hell and Lya, our observations 
would support a an intrinsic Hell/Lya line ratio much closer 
to ^ 0.05 — 0.1, thus strongly favouring a PopIII-like popula¬ 
tion. 

A key question, of course, is whether it is even possible 
or expected to observe Lyman-a coming from PopIII stars 
alone, even if such line is ultra-luminous. The most massive 
PopIII stars should be short-lived (a few Myr), and, without 
any previous contribution to ionise their surroundings from 
e.g. neighbour star clusters or other nearby proto-galaxies, 
the most massive PopIII stars would have to be able to emit 
enough ionising photons to produce a n ionised sphere large r 
than 1 Mpc after less than a few Myrs ( |Cen & Haiman|2000| ), 
before the most massive stars reach the supernovae phase 
and likely start enriching the local environment. However, 
such process (for a single PopIII population in full isolation, 
and fully surrounded by neutral Hydrogen) should take at 
least ^ 3 Myrs to happen: this is simply set by the speed 
of light. However, if neighbouring sources (either PopIII or 
PopII stars) have already contributed towards ionising a local 
bubble, and if PopIII star formation can proceed in a wave-like 
pattern, likely from the highest density regions to the lowest 
densities, by the time later PopIII stars form (still in pristine 
gas which was not contaminated due to being sufficiently far 
away), they will be in ideal conditions to be directly observed 
in Lya. Thus, it may be much more likely to observe po¬ 
tentially composite populations than to observe pure PopIII 
stellar populations. Furthermore, and despite the nature of the 
stellar populations, it is likely that the observability of very 
luminous Lya emitters is strongly favoured in complex sys¬ 
tems which already have older stellar populations like CR7 
(that were able to ionise local bubbles before and thus allow¬ 
ing for strong Lya emission from young stellar populations to 
be observable). 

The detection of the signatures from the first generation of 
stars would provide the first glimpse into the general proper¬ 
ties of these stars, such as their initial mass function (IMF). 
The mass of stars is a fundamental property, as it determines 
the evolution, the chemical enrichment of their surroundings 
and their faith. Theoretical work predicts that th e first gen¬ 
eration of sta rs are very massive, up to 1000 Mq ( [Brornm & 


Larson||2004| ). More recent work is aiming to begin to sam- 


ple the potential IMF of PopIII stars, finding a distribution 
which is very fia t or top-heavy, with masses ranging from 10 
to 1000 Mq (e.g. |Hirano et al.^OlA ), although the PopIII IMF 
is extremely uncertain (e.g. |Greif et al.|20TT] ) and only obser¬ 
vations will be able to constrain it. 


7. CONCLUSIONS 

We presented the spectroscopic follow-up of the two most 
luminous z ^ 6.6 Lya candidates in the COSMOS field 
(LLya - 3 - 9 X 10^^ erg s"^) : ‘MASOSA and ‘C R7’. 
These sources were identified in [Matthee et al.| ( |2015| ), re¬ 
vealing that such luminous sources are much more com¬ 
mon than previously thought and have number densities of 
^ 1.5 X 10“^Mpc“^. Our main results are: 


VLT, and DEIMOS on Keck, to confirm both candi¬ 
dates beyond any doubt. We find redshifts of 2 ; = 
6.604 and z = 6.541 for ‘CR7’ and ‘MASOSA, re¬ 
spectively. ‘CR7’ has an observed Lya luminosity 
of io^3-^3=tO-0^ergs“^ (^ 3x more luminous than 
Himiko) and is the most luminous Lya emitter ever 
found at the epoch of re-ionisation. 

• MAS OS A has a strong detection in Lya, with very 
high Lya EW (EWq > 200 A), implying very low stel¬ 
lar mass and a likely extreme stellar population. It is, 
nonetheless, undetected in all other available bands and 
Lya is also rather compact. 

• CR7, with a narrow Lya line with 266±15kms“^ 
FWHM, is detected in the NIR (rest-frame UV), with 
/3 = —2.3±0.08, an excess in J band, and it is strongly 
detected in IRAC/Spitzer. 

• We detect a strong Heiil640A narrow emission line 
with both X-SHOOTER and SINFONI in CR7 (imply¬ 
ing 2 ; = 6.600), which is enough to explain the clear 
excess seen in the J band photometry. We find no 
other emission lines from the UV to the NIR in our 
X-SHOOTER spectra. No AGN line is seen, nor any 
signatures of WR stars, as the Helll640A emission line 
is narrow (EWHM = 130 ± 30 km s-^). The Heiil640A 
emission line implies that we are seeing the peak of Lya 
emission redshifted by ±160kms“^ and thus that we 
are only seeing the red wing of the Lya (the intrinsic 
Lya fiux is thus likely much higher than seen), or we 
are witnessing an outflow. 

• The AGN and WR stars interpretation of the nature of 
CR7 are strongly disfavoured. An alternative interpre¬ 
tation is that the source hosts a direct collapse black 
hole, although the lack of broad emission lines and the 
lack of X-ray detection also disfavours this interpreta¬ 
tion. Given all the current data, we conclude that CR7 
may host an unseen, extreme stellar population and it is 
therefore the strongest candidate for PopIII-like stellar 
population found so far. 

• We find that CR7 cannot be described only by a PopIII 
stellar population, particularly due to the very strong 
IRAC detections. Our best interpretation of the full data 
(spectroscopy and photometry), which is fully consis¬ 
tent with many theoretical predictions, is a combination 
of a PopIII-like stellar population, which dominates the 
rest-frame UV and the emission lines, and an older, 
likely metal enriched stellar population, which is red, 
and that completely dominates the mass of the system. 
This interpretation fits remarkably well with high res¬ 
olution i75TAVPC3 imaging that reveals two red com¬ 
ponents, each about 5 kpc away from the peak of the 
rest-frame UV, Lya and Hell 1640 A emission and that 
have the fluxes predicted by our SED fitting. 


We may be witnessing, for the first time, direct evidence for 
the occurrence of waves of PopIII-like star formation which 
could happen from an original star cluster outwards (result¬ 
ing from strong feedback w hich can delay PopIII st ar forma¬ 


tion), as suggested by e.g. Tornatore et al, 


scenario, the reddest clump m CR7 (C, see 
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which formed first (reddest and oldest), was likely responsi¬ 
ble for not only starting to ionise a local bubble, but also for 
producing copious amounts of Lyman-Werner radiation and 
feedback that may have prevented star-formation to occur in 
the vicinity of clump C. Star-formation likely proceeded to the 
second clump once stellar feedback from C declined (B, see 
Figure [7|), with similar effects (preventing star-formation out¬ 
side su3i region, but further ionising a local bubble), and we 
are observing source A (see Figure [ 7 ]) at the right time to see 
an intense PopIII-like star formation episode. Most impor¬ 
tantly, this scenario also provides a very simple explanation 
of why hya photons can easily escape the CR7 galaxy, as we 
see a significant amount of older stars which were able to emit 
a significant amount of UV photons for a few hundred mil¬ 
lion years before observations, enough to ionise a bubble of 
> 1 Mpc around the source. Given the strong Lyman-Werner 
fiux that we can infer is coming from the bright UV core, 
it is very likely that any previous episodes of star-formation 
could have prevented the gas around those star-forming re¬ 
gions to form stars. Furthermore, we note that while radia¬ 
tion is able to affect the surroundings as it travels fast, signifi¬ 
cant metal en richment is very inefficient on scales larger than 
^ Ikpc (e.g. IScannapieco et al.|2003[|Tomatore et al.|2007[ 

Ritter et al.|2014| ), both due to the larger time-scales for metals 

(from supernovae) to travel outwards (compared to the speed 
of light), but also due to the continued infall of metal-free 
gas from the cosmic web. It is therefore likely that in some 
cases, scales beyond l-2kpc of previous star formation activ¬ 
ity can easil y have the pristine gas ne cessary to allow PopIII 
to form (e.g. [Ritter et al.|20T2||2Q14| ) even at ^ ^ 6.6, and to 
be detectabk at redshifts ev^be l ow ^ r\j 5 (e.g. I cannapiec'ol 
|et al.|2003||Tornatore et al.|2Q07| ). |Tornatore et al.| ( |2007| ), lor 
example, predict that the peak of PopIII star formation rate 
density to occur at 2 ; ^ 6 — 8. 

The spectroscopic confirmation of MASOSA and CR7, 
along with the high S/N spectroscopic and photometric data, 
allowed us to have a first glimpse into sources likely similar to 
Himiko and brighter, that are much more common than previ¬ 
ously exp ected and have a remar kable nature. The follow-up 
of the full [Matthee et al.| ( |2015| ) Lya sources at even higher 
luminosities found in the SA22 field will allow us to ex¬ 
plore even more the diversity and nature of such unique tar¬ 
gets. Such luminous Lya emitters are the ideal first targets 
for JWST, particularly due to the likely very high EW and 
bright optical rest-frame emission lines, which may not be 
restricted to bright [Olll]5007, H/3 and Ha, but actually in¬ 
clude Heii4686, Hei4471, Hei5016 and Hei5876. In the case 
of CR7, a composite of a PopIII-like stellar population and 
a likely enriched stellar population which is physically sep¬ 
arated by ^ 5 kpc (Figure [ 7 ]) is currently strongly favoured, 
and it is possible that similar stellar populations will be found 
in the other Lya emitters. We may have found an actual popu¬ 
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lation with clear signatures of PopIII-like stars, besides CR7. 
JWST will, in only a short exposure time, clearly show if the 
rest-frame spectra is made up of only He-i-H lines, confirm¬ 
ing PopIII beyond any doubt, or if [Olll] is also present and 
exactly where each of the lines is coming from. 
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